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Abstract Results of a search for new phenomena in final states with an energetic jet and large missing 
transverse momentum are reported. The search uses 20.3 fb _1 of y/s = 8 TeV data collected in 2012 with 
the ATLAS detector at the LHC. Events are required to have at least one jet with px > 120 GeV and 
no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements 
between E™ 1SS > 150 GeV and E™ 1SS > 700 GeV. Good agreement is observed between the number of 
events in data and Standard Model expectations. The results are translated into exclusion limits on models 
with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, 
or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on 
the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are 
presented. 


1 Introduction 

Events with an energetic jet and large missing transverse momentum in the final state constitute a clean 
and distinctive signature in searches for new physics beyond the Standard Model (SM) at colliders. Such 
signatures are referred to as monojet-like in this paper. In particular, monojet-like (as well as monophoton 
and mono-IT/^) final states have been studied [I— 20] in the context of searches for supersymmetry (SUSY), 
large extra spatial dimensions (LED), and the search for weakly interacting massive particles (WIMPs) as 
candidates for dark matter (DM). 

The Arkani-Hamed, Dimopoulos, and Dvali (ADD) model for LED [21] explains the large difference 
between the electroweak unification scale at 0( 10 1 2 * 4 ) GeV and the Planck scale M P1 ~ O(10 19 ) GeV by 
postulating the presence of n extra spatial dimensions of size R , and defining a fundamental Planck scale in 

4 + n dimensions, Mp, given by Mpi 2 ~ Mp 2+n i? rl . An appropriate choice of R for a given n yields a value 
of Mj-j at the electroweak scale. The extra spatial dimensions are compactified, resulting in a Kaluza-Klein 
tower of massive graviton modes. If produced in high-energy collisions in association with an energetic jet, 
these graviton modes escape detection leading to a monojet-like signature in the final state. 

A non-baryonic DM component in the universe is commonly used to explain a range of astrophysical 
measurements (see, for example, Ref. [22] for a review). Since none of the known SM particles are adequate 
DM candidates, the existence of a new particle is often hypothesized. Weakly interacting massive particles 
are one such class of particle candidates that can be searched for at the LHC [23]. They are expected to 
couple to SM particles through a generic weak interaction, which could be the weak interaction of the SM 
or a new type of interaction. Such a new particle would result in the correct relic density values for non- 
relativistic matter in the early universe [24], as measured by the PLANCK [25] and WMAP [26] satellites, 
if its mass is between a few GeV and a TeV and if it has electroweak-scale interaction cross sections. Many 
new particle physics models such as SUSY [27—35] also predict WIMPs. 

Because WIMPs interact so weakly that they do not deposit energy in the calorimeter, their production 
leads to signatures with missing transverse momentum (pt? 1188 ), the magnitude of which is called E™ 1SS . 
Here, WIMPs are assumed to be produced in pairs, and the events are identified via the presence of an 
energetic jet from initial-state radiation (ISR) [36-39] yielding large E™ 1SS . 

The interaction of WIMPs with SM particles is described as a contact interaction using an effective 
field theory (EFT) approach, mediated by a single new heavy particle or particles with mass too large to 
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Fig. 1 Feynman diagrams for the production of weakly interacting massive particle pairs XX associated with a jet from 
initial-state radiation of a gluon, g. (a) A contact interaction described with effective operators, (b) A simplified model with 
a Z' boson. 


be produced directly at the LHC (see Fig. 1(a)). It is assumed here that the DM particle is either a Dirac 


Table 1 Effective interactions coupling WIMPs to Standard Model quarks or gluons, following the formalism in Ref. [40], 
where M* is the suppression scale of the interaction. Operators starting with a D describe Dirac fermion WIMPs, the ones 
starting with a C are for scalar WIMPs and G“„ is the colour field-strength tensor. 


Name 

Initial state 

Type 

Operator 

Cl 

<?<? 

scalar 

m| X r X99 

C5 

99 

scalar 

(G“„) 2 

Dl 

99 

scalar 

Tjkxxqq 

D5 

99 

vector 


D8 

99 

axial-vector 

1^X7 m 7 5 X97m7 5 9 

D9 

99 

tensor 

■j^jx^xq^^q 

Dll 

99 

scalar 

4j%XXas(G“J 2 


fermion or a scalar x; the only difference for Majorana fermions is that certain interactions are not allowed 
and that the cross sections for the allowed interactions are larger by a factor of four. Seven interactions are 
considered (see Table 1), namely those described by the operators Cl, C5, Dl, D5, D8, D9, Dll, following 
the naming scheme in Ref. [40]. These operators describe different bilinear quark couplings to WIMPs, 
qq —> XX) except for C5 and Dll, which describe the coupling to gluons, gg —> \X- The operators for 
Dirac fermions and scalars in Ref. [40] fall into six categories with characteristic E™ 1SS spectral shapes. The 
representative set of operators for these six categories are Cl, C5, Dl, D5, D9, and Dll, while D8 falls 
into the same category as D5 but is listed explicitly in Table 1 because it is often used to convert LHC 
results into limits on DM pair production. In the operator definitions in Table 1 , M* is the suppression scale 
of the interaction, after integrating out the heavy mediator particles. The use of a contact interaction to 
produce WIMP pairs via heavy mediators is considered conservative because it rarely overestimates cross 
sections when applied to a specific scenario for physics beyond the SM. Cases where this approach is indeed 
optimistic are studied in Refs. [39, I 15] . Despite the caveats related to the validity of the EFT approach 
(see Appendix A), this formalism is used here, as it provides a framework for comparing LHC results to 
existing direct or indirect DM searches. Within this framework, interactions of SM and DM particles are 
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described by only two parameters, the suppression scale M* and the DM particle mass m x . Besides the EFT 
operators, the pair production of WIMPs is also investigated within a so-called simplified model, where a 
pair of WIMPs couples to a pair of quarks explicitly via a new mediator particle, a new vector boson Z' 
(see Fig. 1(b)). 

In gauge-mediated SUSY-breaking (GMSB) scenarios [46-51], the gravitino G (spin-3/2 superpartner 
of the graviton) is often the lightest supersymmetric particle and a potential candidate for DM. Its mass 
is related to the SUSY-breaking scale \/~F and M P1 via oc F/M P \ [52]. At hadron colliders, in low-scale 
SUSY-breaking scenarios with very light gravitinos, the cross section for associated production of gravitino- 
squark (pp —¥ Gq + X ) and gravitino-gluino (pp —¥ Gg + X ) processes are relatively large [53], since the 
cross section depends on rriQ as a ~ 1/m^,. The decay of the gluino or squark into a gravitino and a gluon 
(g —> Gg) or a gravitino and a quark (q —> Gq), respectively, dominates [53]. The final state is characterized 
by the presence of a pair of gravitinos that escape detection and an energetic jet, leading to a monojet-like 
topology. Previous studies at colliders [15,54] considered the production of gravitinos in association with 
a photon or a jet and assumed extremely heavy squarks and gluinos. Within this approximation, a lower 
limit for the gravitino mass of m,Q > 1.37 x 10 -5 eV was established. 

The study of the properties of the Higgs boson discovered by the ATLAS and CMS experiments [55,56] 
does not exclude a sizeable branching ratio for its decay to invisible particles. It also opens up the question 
of whether a Higgs-like scalar field plays an important role in describing the interaction between dark and 
ordinary matter in the universe. In particular, a sizeable branching ratio to invisible particles could be 
interpreted in terms of the production of DM. Results from LEP [57] excluded an invisibly decaying Higgs 
boson, produced in association with a Z boson, for a boson mass (m#) below 114.4 GeV. The strongest 
direct bounds from the LHC experiments on the branching ratio for the Higgs invisible decay mode [58,59] 
set upper limits of 58%-65% at 95% confidence level (CL), based on the final state in which the Higgs 
boson is produced either in association with a Z boson or via vector-boson fusion processes. In this analysis, 
the monojet-like final state is used to search for the production of an invisibly decaying boson with SM 
Higgs-like properties and a mass in the range between 115 GeV and 300 GeV. 

The paper is organized as follows. The ATLAS detector is described in the next section. Section 3 provides 
details of the simulations used in the analysis for background and signal processes. Section 4 discusses the 
reconstruction of jets, leptons and E™ 1SS , while Sect. 5 describes the event selection. The estimation of 
background contributions and the study of systematic uncertainties are discussed in Sects. 6 and 7. The 
results are presented in Sect. 8, and are interpreted in terms of the search for ADD LED, WIMP pair 
production, the production of very light gravitinos in GMSB scenarios, and the production of an invisibly 
decaying Higgs-like boson. Finally, Sect. 9 is devoted to the conclusions. 


2 Experimental setup 

The ATLAS detector [60] covers almost the whole solid angle 1 around the collision point with layers of 
tracking detectors, calorimeters and muon chambers. The ATLAS inner detector (ID) has full coverage in 
4> and covers the pseudorapidity range \r/\ < 2.5. It consists of a silicon pixel detector, a silicon microstrip 
detector, and a straw tube tracker which also measures transition radiation for particle identification, all 
immersed in a 2 T axial magnetic field produced by a solenoid. 

High-granularity liquid-argon (LAr) electromagnetic sampling calorimeters, with excellent energy and 
position resolution, cover the pseudorapidity range |r?| < 3.2. The hadronic calorimetry in the range |r 7 | <1.7 
is provided by a scintillator-tile calorimeter, consisting of a large barrel and two smaller extended barrel 
cylinders, one on either side of the central barrel. In the endcaps (I?/ > 1.5), LAr hadronic calorimeters 
match the outer |r;| limits of the endcap electromagnetic calorimeters. The LAr forward calorimeters provide 
both the electromagnetic and hadronic energy measurements, and extend the coverage to |r?| < 4.9. 

The muon spectrometer measures the deflection of muons in the magnetic field provided by large su¬ 
perconducting air-core toroid magnets in the pseudorapidity range I?/ < 2.7, instrumented with separate 
trigger and high-precision tracking chambers. Over most of the g range, a measurement of the track coor¬ 
dinates in the principal bending direction of the magnetic field is provided by monitored drift tubes. At 

1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the 
detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis 
points upward. The azimuthal angle <f) is measured around the beam axis, and the polar angle 6 is measured with respect 
to the z-axis. We define transverse energy as /V[- = E sin0, transverse momentum as fry = p si rid, and pseudorapidity as 
g = —ln[tan(0/2)[. 
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large pseudorapidities, cathode strip chambers with higher granularity are used in the innermost plane over 
2.0 < |r/| < 2.7. The muon trigger system covers the pseudorapidity range \r/\ < 2.4. 

The data are collected using an online three-level trigger system [61] that selects events of interest and 
reduces the event rate from several MHz to about 400 Hz for recording and offline processing. 


3 Monte Carlo simulation 

Monte Carlo (MC) simulated event samples are used to compute detector acceptance and reconstruction 
efficiencies, determine signal and background contributions, and estimate systematic uncertainties on the 
final results. 


3.1 Background simulation 

The expected background to the monojet-like signature is dominated by Z(—> uv)+ jets and W+jets produc¬ 
tion (with W{—> rr'j+jets being the dominant among the VF+jets backgrounds), and includes small contri¬ 
butions from Zj 7 *(—>■ i + l~)+]ets (£ = e,/x,r), multijet, it, single-top, and diboson (WW,WZ, ZZ,W-y,Z"/) 
processes. 

Samples of simulated W+jets and -Z'+jets production events are generated using SHERPA-1.4.1 [62], 
including leading-order (LO) matrix elements for up to five partons in the final state and assuming mas¬ 
sive 6 /c-quarks, with CT10 [63] parton distribution functions (PDF) of the proton. The MC expectations 
are initially normalized to next-to-next-to-leading-order (NNLO) perturbative QCD (pQCD) predictions 
according to DYNNLO [64,65] using MSTW2008 90% CL NNLO PDF sets [ 66 ]. The production of top- 
quark pairs ( tt) is simulated using the MC@NLO-4.06 [67, 68 ] MC generator with parton showers and 
underlying-event modelling as implemented in HERWIG-6.5.20 [69,70] plus JIMMY [71]. Single-top pro¬ 
duction samples are generated with MC@NLO [72] for the s- and Wt -channel [73], while AcerMC-v3.8 [74] 
is used for single-top production in the t-channel. A top-quark mass of 172.5 GeV is used consistently. 
The AUET2C and AUET2B [75] set of optimised parameters for the underlying event description are 
used for tt and single-top processes, which use CT10 and CTEQ6L1 [76] PDF, respectively. Approximate 
NNLO+NNLL (next-to-next-to-leading-logarithm) pQCD cross sections, as determined in TOP++2.0 [77], 
are used in the normalization of the tt [78] and Wt [79] samples. Multijet and 7 +jet samples are generated 
using the PYTHIA-8.165 program [80] with CT10 PDF. Finally, diboson samples {WW , WZ, ZZ, W 7 
and Z 7 production) are generated using SHERPA with CT10 PDF and are normalized to NLO pQCD 
predictions [81]. 


3.2 Signal simulation 

Simulated samples for the ADD LED model with different number of extra dimensions in the range n =2— 
6 and Mjj in the range 2-5 TeV are generated using PYTHIA-8.165 with CT10 PDF. Renormalization 
and factorization scales are set to \j 1/2 x tuq + where me is the graviton mass and pt denotes the 
transverse momentum of the recoiling parton. 

The effective field theory of WIMP pair production is implemented in MADGRAPH5-vl.5.2 [82], taken 
from Ref. [40]. The WIMP pair production plus one or two additional partons from ISR is simulated in 
two ways. For all operators, samples are generated requiring at least one parton with a minimum px of 
80 GeV. Studies simulating up to three additional partons along with the WIMP pair showed no difference 
in kinematic distributions when compared to the samples with up to two additional partons. 

Only initial states of gluons and the four lightest quarks are considered, assuming equal coupling 
strengths for all quark flavours to the WIMPs. The mass of the charm quark is most relevant for the 
cross sections of the operator D1 (see Table 1) and it is set to 1.42 GeV. The generated events are interfaced 
to PYTHIA-6.426 [83] for parton showering and hadronization. The MLM prescription [84] is used for 
matching the matrix-element calculations of MADGRAPH5 to the parton shower evolution of PYTHIA- 6 . 
The samples are subsequently reweighted to the MSTW2008LO [ 66 ] PDF set using LHAPDF [85]. The 
MADGRAPH5 default choice for the renormalization and factorization scales is used. The scales are set to 
the geometric average of m 2 + p\ for the two WIMPs, where m is the mass of the particles. Events with 
WIMP masses between 10 GeV and 1300 GeV are simulated for six different effective operators (Cl, C5, 
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Dl, D5, D9, Dll). The WIMPs are taken to be either Dirac fermions (D operators) or scalars (C operators), 
and the pair-production cross section is calculated at LO. To study the transition between the effective field 
theory and a physical renormalizable model for Dirac fermion WIMPs coupling to Standard Model particles 
via a new mediator particle Z', a simplified model is generated in MADGRAPH5. For each WIMP mass 
point, mediator particle masses M med between 50 GeV and 30 TeV are considered, each for two values of 
the mediator particle width (F = M me d/3 and M me ^/8n). 

Simulated samples for gravitino production in association with a gluino or a squark in the final state, 
pp —> Gg + X and pp — > Gq + X , are generated using LO matrix elements in MADGRAPH4.4 [86] interfaced 
with PYTHIA-6.426 and using CTEQ6L1 PDF. The narrow-width approximation for the gluino and squark 
decays g — > gG and q —> qG is assumed. The renormalization and factorization scales are set to the average 
of the mass of the final-state particles involved in the hard interaction (mg + m^/g )/2 ~ rrig/g/2. Values for 
mg in the range between 10 “ 3 eV and 10 -5 eV are considered for squark and gluino masses in the range 
50 GeV to 2.6 TeV. 

Finally, MC simulated samples for the production of a Higgs boson are generated including the gg —> H, 
VV —¥ H (V = W, Z), and VH production channels. Masses for the boson in the range between 115 GeV 
and 300 GeV are considered. This Higgs boson is assumed to be produced as predicted in the Standard 
Model but unlike the SM Higgs it may decay into invisible particles at a significant rate. The signal is 
modelled using POWHEG-r2262 [87-89], which calculates separately the gg —> H and VV —> H production 
mechanisms with NLO pQCD matrix elements. The description of the Higgs boson px spectrum in the 
gg —> H process follows the calculation in Ref. [90], which includes NLO + NNLL corrections. The effects 
of finite quark masses are also taken into account [91]. For gg —» H and VV —> H processes, POWHEG 
is interfaced to PYTHIA-8.165 for showering and hadronization. For ZH and WH processes, POWHEG 
interfaced to HERWIG++ [92] is used and the Z/W bosons are forced to decay to a pair of quarks. The 
invisible decay of the Higgs-like boson is simulated by forcing the boson to decay to two Z bosons, which 
are then forced to decay to neutrinos. Signal samples are generated with renormalization and factorization 
scales set to \J (m#) 2 + (p^f) 2 . The Higgs boson production cross sections, as well as their uncertainties, 
are taken from Refs. [93,94]. For the gg —» H process, cross-section calculations at NNLO+NNLL accuracy 
[95-98] in pQCD are used and NLO electroweak corrections [99,100] are included. The cross sections for 
VV —> H processes are calculated with full NLO pQCD and electroweak corrections [101-103]. The cross 
sections for the associated production (WH and ZH) are calculated at NNLO [104] in pQCD, and include 
NLO electroweak corrections [105]. 

Differing pileup (multiple proton—proton interactions in the same or neighbouring bunch-crossings) 
conditions as a function of the instantaneous luminosity are taken into account by overlaying simulated 
minimum-bias events generated with PYTHIA-8 onto the hard-scattering process. The MC-generated sam¬ 
ples are processed either with a full ATLAS detector simulation [106] based on the GEANT4 program [107] 
or a fast simulation of the response of the electromagnetic and hadronic calorimeters [108] and of the trigger 
system. The results based on fast simulation are validated against fully simulated samples and the difference 
is found to be negligible. The simulated events are reconstructed and analysed with the same analysis chain 
as for the data, using the same trigger and event selection criteria. 


4 Reconstruction of physics objects 

Jets are defined using the anti-fc* jet algorithm [109] with the radius parameter R = 0.4. Energy depositions 
reconstructed as clusters in the calorimeter are the inputs to the jet algorithm. The measured jet px is 
corrected for detector effects, including the non-compensating character of the calorimeter, by weighting 
energy deposits arising from electromagnetic and hadronic showers differently. In addition, jets are corrected 
for contributions from pileup, as described in Ref. [110]. Jets with corrected px > 30 GeV and |p| < 4.5 
are considered in the analysis. Jets with |p| < 2.5 containing a b-hadron are identified using a neural-net- 
based algorithm [111] with an efficiency of 80% and a rejection factor of 30 (3) against jets originating 
from fragmentation of light quarks or gluons (jets containing a c-hadron), as determined using tt simulated 
events. 

The presence of leptons (muons or electrons) in the final state is used in the analysis to define control 
samples and to reject background contributions in the signal regions (see Sects. 5 and 6). Muon candidates 
are formed by combining information from the muon spectrometer and inner tracking detectors as described 
in Ref. [112] and are required to have px > 7 GeV and \r/\ < 2.5. In addition, muons are required to be 
isolated: the sum of the transverse momenta of the tracks not associated with the muon in a cone of size 
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AR = ( Ar)) 2 + ( A<f> ) 2 = 0.2 around the muon direction is required to be less than 1.8 GeV. The muon 
Px requirement is increased to px > 20 GeV to define the JE(—»/m)+ jets and Z/ 7 *(—>■ p + fj,~) + jets control 
regions. 

Electron candidates are initially required to have px > 7 GeV and |p| < 2.47, and to pass the medium 
electron shower shape and track selection criteria described in Ref. [113], which are reoptimized for 2012 
data. Overlaps between identified electrons and jets in the final state are resolved. Jets are discarded if 
their separation AR from an identified electron is less than 0.2. The electron px requirement is increased 
to px > 20 GeV and the transition region between calorimeter sections 1.37 < |p| < 1.52 is excluded to 
reconstruct Z and W boson candidates in the Z/q*(—e + e - )+jets and W(—r e^)+jets control regions, 
respectively. The electron requirements are further tightened for the W{—> ez'j+jets control sample to 
constrain the irreducible Z(—> vv)+ jets background contribution (see below). In this case, electrons are 
selected to pass tight [113] electron shower shape and track selection criteria, their px threshold is raised 
to 25 GeV, and they are required to be isolated: the sum of the transverse momenta of the tracks not 
associated with the electron in a cone of radius AR = 0.3 around the electron direction is required to be 
less than 5% of the electron px. An identical isolation criterion, based on the calorimeter energy deposits 
not associated with the electron, is also applied. 

The E^ 1SS is reconstructed using all energy deposits in the calorimeter up to pseudorapidity |p| = 4.9. 
Clusters associated with either electrons or photons with px > 10 GeV and those associated with jets with 
Px > 20 GeV make use of the corresponding calibrations for these objects. Softer jets and clusters not 
associated with these objects are calibrated using both calorimeter and tracking information [114]. 


5 Event selection 

The data sample considered in this paper corresponds to a total integrated luminosity of 20.3 fb _1 . The 
uncertainty on the integrated luminosity is 2.8%, as estimated following the same methodology as detailed in 
Ref. [115]. The data were selected online using a trigger logic that selects events with E™ 1SS above 80 GeV, 
as computed at the final stage of the three-level trigger system [61]. With respect to the final analysis 
requirements, the trigger selection is fully efficient for E™ 1SS > 150 GeV, as determined using a data sample 
with muons in the final state. Table 2 summarizes the different event selection criteria applied in the signal 
regions. The following preselection criteria are applied. 

— Events are required to have a reconstructed primary vertex for the interaction consistent with the 
beamspot envelope and to have at least two associated tracks with px > 0.4 GeV; when more than one 
such vertex is found, the vertex with the largest summed Px of the associated tracks is chosen. 

— Events are required to have E™' ss > 150 GeV and at least one jet with px > 30 GeV and |p| < 4.5 in 
the final state. 

— The analysis selects events with a leading jet with px > 120 GeV and |p| < 2.0. Monojet-like topologies 
in the final state are selected by requiring the leading-jet px and the E™ 1SB to satisfy pt/E™ 1ss > 0.5. 
An additional requirement on the azimuthal separation A^>(jet, Pt?" ss ) > 1.0 between the direction of 
the missing transverse momentum and that of each of the selected jets is imposed. This requirement 
reduces the multijet background contribution where the large E™ 1SS originates mainly from jet energy 
mismeasurement. 

— Events are rejected if they contain any jet with px > 20 GeV and |p| < 4.5 that presents an electro¬ 
magnetic fraction in the calorimeter, calorimeter sampling fraction, or charged fraction 2 (for jets with 
|p| < 2.5) inconsistent with the requirement that they originate from a proton-proton collision [116]. In 
the case of the leading (highest px) jet in the event, the requirements are tightened to reject remaining 
contributions from beam-related backgrounds and cosmic rays. Events are also rejected if any of the jets 
is reconstructed close to known partially instrumented regions of the calorimeter. Additional require¬ 
ments based on the timing and the pulse shape of the cells in the calorimeter are applied to suppress 
coherent noise and electronic noise bursts in the calorimeter producing anomalous energy deposits [117]; 
these requirements have a negligible effect on the signal efficiency. 

— Events with muons or electrons with px > 7 GeV are vetoed. In addition, events with isolated tracks 
with px > 10 GeV and |p| < 2.5 are vetoed to reduce background from non-identified leptons (e, p or 

2 The charged fraction is defined as / c h = V ^trackjet w } U!re ^ptrack,jet ^ t } |e sca j ar sum Q f the transverse momenta 

of tracks associated with the primary vertex within a cone of radius R = 0.4 around the jet axis, and is the transverse 
momentum as determined from calorimetric measurements. 
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r) in the final state. The track isolation is defined such that there must be no additional track with 
Pt > 3 GeV within a cone of radius 0.4 around it. 

Different signal regions (SR1-SR9) are considered with increasing E™ 1SS thresholds from 150 GeV to 
700 GeV. 


Table 2 Event selection criteria applied for the selection of monojet-like signal regions, SR1-SR9. 


Selection criteria 


Preselection 


Primary vertex 

E“ iss > 150 GeV 

Jet quality requirements 

At least one jet with px > 30 GeV and \rj\ < 4.5 
Lepton and isolated track vetoes 

Monojet-like selection 

The leading jet with px > 120 GeV and ?/ < 2.0 
Leading jet px/-E™ lss > 0.5 
A0(jet,p T miss ) > 1.0 

Signal region SRI SR2 SR3 SR4 SR5 SR6 SR7 SR8 SR9 

Minimum Kj? iss [GeV] 150 200 250 300 350 400 500 600 700 


6 Background estimation 

The W -t-jets and Z(—t vv)+jets backgrounds are estimated using MC event samples normalized using data in 
selected control regions. In particular, the dominant Z(—t vv)+jets background contribution is constrained 
using a combination of estimates from IV+jets and Z +jets control regions. The remaining SM backgrounds 
from Z/ 7 *(—t t + £ _ )+jets, tt, single top, and dibosons are determined using MC simulated samples, while 
the multijet background contribution is extracted from data. In the case of the tt background process, 
which contributes to both the signal and IV+jets control regions, dedicated control samples are defined to 
validate the MC normalization and to estimate systematic uncertainties. Finally, the potential contributions 
from beam-related background and cosmic rays are estimated in data using jet timing information. The 
methodology and the samples used for estimating the background are summarised in Table 3. The details 
are given in the following sections. 


Table 3 Summary of the methods and control samples used to constrain the different background contributions in the 
signal regions. 


Background process 

Method 

Control sample 

Z(—>• i/p)+jets 

MC and control samples in data 

Z/ 7*(-> t+t~), W(-»- tv) {t = e,p) 

W(—t ev)+jets 

MC and control samples in data 

W (— > ev) (loose) 

W (—> Ti/)+jets 

MC and control samples in data 

W (— > ev) (loose) 

W(—> /nd+jets 

MC and control samples in data 

W(-+ fiv) 

Z/ 7 *(—>■ £ + t - )+jets (t = e,p,r) 

MC-only 


tt, single top 

MC-only 


Diboson 

MC-only 


Multijets 

data-driven 


Non-collision 

data-driven 

















6.1 W/Z+jets background 


In the analysis, control samples in data, with identified electrons or muons in the final state and with 
identical requirements on the jet px and E™ 1SS , are used to determine the W(—> £v)+)ets (£ = e,/r, t) 
and Z(—> j/i/)+jets electroweak background contributions. This reduces significantly the relatively large 
theoretical and experimental systematic uncertainties, of the order of 20%-40%, associated with purely MC- 
based expectations. The i?x ISS -based online trigger used in the analysis does not include muon information 
in the E™ lss calculation. This allows the collection of W(—> pv)+ jets and Z/^*{—> p + p~)+ jets control 
samples with the same trigger as for the signal regions. This is not the case for the W{—> eiz)+jets and 
Zh*{~> e + e _ )+jets control samples used to constrain the Z{—> vv)+]eis background (see below). 

A VF (—> /r^)+jets control sample is defined using events with a muon withpx > 20 GeV and W transverse 
mass in the range 40 GeV < mx < 100 GeV. The transverse mass mx is defined by the lepton (£) and 

neutrino ( v) px and direction as mp = \j2p^ r pY r {\ — coswhere the ( x,y ) components of the 
neutrino momentum are taken to be the same as the corresponding pr,? 115 " 3 components. Similarly, a Z/^*{—> 
p + p - )+jets control sample is selected, requiring the presence of two muons withpx > 20 GeV and invariant 
mass in the range 66 GeV < < 116 GeV. In the W(—¥ pi^)+jets and Z /^*{—> p + p“)+jets control 

regions, the E™ lss is not corrected for the presence of the muons in the final state, which are considered 
invisible, motivated by the fact that these control regions are used to estimate the irreducible Z(—> i'Pj+jets 
background in the signal regions. 

The W{—> ei/j+jets and Zj 7 *(—» e + e _ )+jets control samples used to constrain the Z{—¥ vv)+]ets back¬ 
ground in the signal regions are collected using online triggers that select events with an electron in the 
final state. The E™ 1SS is corrected by removing the contributions from the electron energy clusters in the 
calorimeters. In the Z/^*{—> e + e _ )4-jets control sample, events are selected with exactly two electrons with 
Px > 20 GeV and dilepton invariant mass in the range 66 GeV < m ee < 116 GeV. In the W{—> e^)+jets con¬ 
trol sample a tight selection is applied: events are selected to have only a single electron with px > 25 GeV, 
transverse mass in the range 40 GeV < mx < 100 GeV, and uncorrected E™ 1SS > 25 GeV. The latter 
requirements suppress background contamination from multijet processes where jets are misidentified as 
electrons. 

A separate W(—t ei/)+jets control sample, collected with the i?x iss -based trigger and looser requirements 
that increase the number of events, is defined to constrain the W(—> e^j+jets and W(—> ri/)-|-jets background 
contributions. In this case, the electron px requirement is reduced to px > 20 GeV and no further cuts on 
electron isolation and mx are applied. In addition, the E™ 1SS calculation in this case is not corrected for the 
presence of the electron or r leptons in the final state, as they contribute to the calorimeter-based E™ 1SS 
calculation in the signal regions. 

Figures 2-5 show, for the SRI monojet-like kinematic selection, some distributions in data in the different 
W+jets and Z+jets control regions compared to MC expectations. In this case, the MC expectations are 
globally normalized to the data in the control regions, using normalization factors as explained below, so 
that a comparison of the shape of the different distributions in data and MC simulation can be made. The 
MC expectations provide a fair description of the shapes in data but present harder E™ lss and leading-jet px 
spectra. This is mainly attributed to an inadequate modelling of the boson px distribution in the IV/Z+jets 
MC samples. 

The data in the control regions and MC-based correction factors, determined from the SHERPA simu¬ 
lation, are used for each of the signal selections (SR1-SR9) to estimate the electroweak background contri¬ 
butions from IT+jets and Z{—¥ vv)+jets processes. As an example, the W(—> pv)+)ets and Z(—t I'i'j+jets 
background contributions to a given signal region, and respectively, are determined 

using the W(—> /xi/)+jets control sample in data according to 
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Fig. 2 Distributions of the measured (a) transverse mass of the identified muon and the missing transverse momentum, 
(b) £vp lss , (c) leading jet p t and (d) jet multiplicity distributions in the piz)-\-]ets control region for the inclusive SRI 

selection, compared to the background expectations. The latter include the global normalization factors extracted from the 
data. Where appropriate, the last bin of the distribution includes overflows. The lower panels represent the ratio of data to 
MC expectations. The error bands in the ratios include the statistical and experimental uncertainties on the background 
expectations. 


control denote, in the control region, the number of W (—»■ /xi/)+jets candidates in data and W/Z -\-jets 

MC simulation, and the non -W/Z background contribution, respectively. The ^w(-^pl) Z c 0 ntTo\ term refers 
mainly to top-quark and diboson processes, but also includes contributions from multijet processes de¬ 
termined using data. Finally, £j eto , and £trg account for possible data-MC differences in the lepton 
identification, lepton vetoes, and trigger efficiencies, respectively; they typically depart from unity by less 
than 1%. The MC-to-data normalization factors (the (N^^ hcontml - ^p^ ontrol )/^f-».^) .control 
term in Eq. (1)) for each process vary between about 0.9 and 0.6 as the required minimum E™ 1SS increases 
from 150 GeV to 700 GeV, and account for the tendency of the MC expectations for W/Z +jets processes 
to exceed the data in the control regions (see, for example, Fig. 2). Similarly, bin-by-bin correction factors 
are used to correct the shape of the different distributions in the signal regions. 

As already mentioned, the different background contributions in the signal regions from W{—> £u)+ jets 
processes (with £ = e, p) are constrained using correction factors obtained from the corresponding control 
regions. In the case of the W(—> TZ/)+jets contributions, the correction factors from the W(—» ei/)+jets con¬ 
trol regions are used. For each of the signal regions, four separate sets of correction factors are considered to 
constrain the dominant Z{—¥ vis)+jets background contribution, following Eq. (2), as determined separately 
using Z/ 7 *(—» £ + £~)+]ets and W(—> £^)-|-jets control samples. The four resulting Z(—t ^i^+jets background 
estimations in each signal region are found to be consistent within uncertainties and are statistically com- 
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Fig. 3 Distributions of the measured (a) dilepton invariant mass, (b) E™ laa , (c) leading jet pry and (d) jet multiplicity 
distributions in the Z/ 7*(—>■ )+jets control region for the inclusive SRI selection, compared to the background 

expectations. The latter include the global normalization factors extracted from the data. Where appropriate, the last bin 
of the distribution includes overflows. The lower panels represent the ratio of data to MC expectations. The error bands in 
the ratios include the statistical and experimental uncertainties on the background expectations. 


bined using the BLUE (Best Linear Unbiased Estimate) [118] method, which takes into account correlations 
of systematic uncertainties. 


6.2 Multijet background 

The multijet background with large E™ 1SS mainly originates from the misreconstruction of the energy of a 
jet in the calorimeter and to a lesser extent from the presence of neutrinos in the final state due to heavy- 
flavour decays. The multijet background is determined from data, using a jet smearing method as described 
in Ref. [119], which relies on the assumption that the E™ BS of multijet events is dominated by fluctuations in 
the detector response to jets measured in the data. For the SRI and SR2 selections, the multijet background 
constitutes about 2% and 0.7% of the total background, respectively, and is below 0.5% for the rest of the 
signal regions with higher E™ 1SS thresholds. 


6.3 Non-collision background 

Detector noise, beam-halo and cosmic muons leading to large energy deposits in the calorimeters represent a 
significant portion of data acquired by E™ 1SS triggers. These non-collision backgrounds resemble the topology 
of monojet-like final states and require a dedicated strategy to suppress them. The selection described in 
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Fig. 4 Distributions of the measured (a) transverse mass of the identified electron and the missing transverse momentum, 
(b) Kj? lss , (c) leading jet p-[- and (d) jet multiplicity distributions in the W(—> evj+jets control region for the inclusive SRI 
selection, compared to the background expectations. The latter include the global normalization factors extracted from the 
data. Where appropriate, the last bin of the distribution includes overflows. The lower panels represent the ratio of data to 
MC expectations. The error bands in the ratios include the statistical and experimental uncertainties on the background 
expectations. 


Sect. 5 is expected to maintain the non-collision background below the percent level. The rate of the fake 
jets due to cosmic muons surviving the selection criteria, as measured in dedicated cosmic datasets, is found 
negligible with respect to the rate of data in the monojet-like signal regions. The major source of the non¬ 
collision backgrounds is thus beam-halo muons. Since jets due to collisions are expected to be in time with 
the bunch crossing, an assumption is made that all events containing a leading jet within the out-of-time 
window are due to beam-induced backgrounds. The characteristic shape of the fake jets due to beam-halo 
muons is extracted from signal-region events identified as beam-induced backgrounds based on the spatial 
alignment of the signals in the calorimeter and the muon system [116]. The level of non-collision background 
in the signal region is extracted as 


n ncb 


/v SR v 

J v -10<i<-5 X 


N 


NCB 


iwNCB ’ 

iv -10<t<-5 


(3) 


where A^io<i <-5 denotes the number of events in the signal region with a leading jet in the range 
— 10 ns < t < —5 ns, A r £ f p 0 ^ ;t< _ 5 is the number of identified beam-induced background events there and 
N ncb represents all identified events in the signal region. The results of this study indicate that the non¬ 
collision background in the different signal regions is negligible. 
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Fig. 5 Distributions of the measured (a) dilepton invariant mass, (b) F™ lss , (c) leading jet pry and (d) jet multiplicity 
distributions in the Z/ 7*(—t e+e — )+jets control region for the inclusive SRI selection, compared to the background expec¬ 
tations. The latter include the global normalization factors extracted from the data. Where appropriate, the last bin of the 
distribution includes overflows. The lower panels represent the ratio of data to MC expectations. The error bands in the 
ratios include the statistical and experimental uncertainties on the background expectations. 


7 Systematic uncertainties 

Several sources of systematic uncertainty are considered in the determination of the background contribu¬ 
tions. Uncertainties on the absolute jet energy scale and resolution [110] translate into an uncertainty on 
the total background which varies from 0.2% for SRI and 1% for SR7 to 3% for SR9. Uncertainties on the 
E™ 1SS reconstruction introduce an uncertainty on the total background which varies from 0.2% for SRI and 
0.7% for SR7 to 1% for SR9. Uncertainties of the order of l%-2% on the simulated lepton identification and 
reconstruction efficiencies, energy/momentum scale and resolution, and a 0.5%-l% uncertainty on the track 
isolation efficiency translate, altogether, into a 1.4%, 1.5%, and 2% uncertainty in the total background 
for the SRI, SR7, and SR9 selections, respectively. Uncertainties of the order of 1% on the E™ 1SS trigger 
simulation at low E™ 1SS and on the efficiency of the lepton triggers used to define the electron and muon 
control samples translate into uncertainties on the total background of about 0.1% for SRI and become 
negligible for the rest of the signal regions. 

The top-quark-related background contributions, as determined from MC simulations (see Sect. 3), are 
validated in dedicated validation regions defined similarly to the W(—t ei/)+jets and W{—> [iv)+jets control 
regions with A(f>(pr^ liss y jet) > 0.5 and by requiring the presence of two ^-tagged jets in the final state 
with jet \ri\ < 2.4. The comparison between data and MC expectations in those validation regions leads 
to uncertainties on the top-quark background yields which increase from 20% for SRI to 100% for SR7 
and SR9. This translates into uncertainties in the total background expectations which vary from 0.7% for 
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SRI and 2.7% for SR7 to 4% for SR9. Similarly, uncertainties on the simulated diboson background yields 
include uncertainties in the MC generators and the modelling of parton showers employed, variations in 
the set of parameters that govern the parton showers and the amount of initial- and final-state soft gluon 
radiation, and uncertainties due to the choice of renormalization and factorization scales and PDF. This 
introduces an uncertainty on the diboson background expectation which increases from 20% for SRI to 30% 
for SR7 and 80% for SR9. This results in an uncertainty on the total background of 0.7%, 2.3%, and 3% 
for the SRI, SR7, and SR9 selections, respectively. 

Uncertainties on the W/Z +jets modelling include: variations of the renormalization, factorization, and 
parton-shower matching scales and PDF in the SHERPA W/Z+ jets background samples; and uncertainties 
on the parton-shower model considered. In addition, the effect of NLO electroweak corrections on the W +jets 
to -Z’+jets ratio is taken into account [120-122], Altogether, this translates into an uncertainty on the total 
background of about 1% for SRI and SR7 and 3% for SR9. 

Uncertainties on the multijet and y+jets background contamination of 100% and 50%, respectively, in 
the W{—> ez'j+jets control region, propagated to the Z{—> vu )-|-jets background determination in the signal 
regions, introduce an additional 1% uncertainly on the total background for the SR9 selection. The uncer¬ 
tainty on the multijet background contamination in the signal regions leads to a 2% and 0.7% uncertainty on 
the total background for the SRI and SR2 selections, respectively. Finally, the impact of uncertainty on the 
total integrated luminosity, which partially cancels in the data-driven determination of the SM background, 
is negligible. 

After including statistical uncertainties on the data and MC expectations in control regions and on 
the MC expectations in the signal regions, the total background in the signal regions is determined with 
uncertainties that vary from 2.7% for SRI and 6.2% for SR7 to 14% for SR9. 


7.1 Signal systematic uncertainties 

Several sources of systematic uncertainty on the predicted signal yields are considered for each of the 
models for new physics. The uncertainties are computed separately for each signal region by varying the 
model parameters (see Sect. 8). 

Experimental uncertainties include: those related to the jet and E™ 1SS reconstruction, energy scales and 
resolutions; those on the proton beam energy, as considered by simulating samples with the lower and upper 
allowed values given in Ref. [123]; a 1% uncertainty on the trigger efficiency, affecting only SRI; and the 
2.8% uncertainty on the integrated luminosity. Other uncertainties related to the track veto or the jet quality 
requirements are negligible (< 1%). 

Uncertainties affecting the signal acceptance times efficiency Axe, related to the generation of the 
signal samples, include: uncertainties on the modelling of the initial- and final-state gluon radiation, as 
determined using simulated samples with modified parton-shower parameters, by factors of two and one 
half, that enhance or suppress the parton radiation; uncertainties due to PDF and variations of the a B (mz) 
value employed, as computed from the envelope of CT10, MRST2008LO and NNPDF21LO error sets; and 
the choice of renormalization/factorization scales, and the parton-shower matching scale settings, varied by 
factors of two and one half. 

In addition, theoretical uncertainties on the predicted cross sections, including PDF and renormaliza¬ 
tion/factorization scale uncertainties, are computed separately for the different models. 


8 Results and interpretation 

The data and the expected SM expectations in the different signal regions are presented in Tables 4 and 5. 
In general, good agreement is observed between the data and the SM expectations. The largest difference 
between the number of events in data and the expectations is observed in the signal region SR9, correspond¬ 
ing to a 1.7cr deviation with a p-value of 0.05, consistent with the background-only hypothesis. Figures 6 
and 7 show several measured distributions in data compared to the SM expectations for SRI, and SR7 and 
SR9, respectively. For illustration purposes, the distributions include the impact of different ADD, WIMP, 
and GMSB SUSY scenarios. 

The agreement between the data and the SM expectations for the total number of events in the different 
signal regions is translated into model-independent 90% and 95% confidence level (CL) upper limits on the 
visible cross section, defined as the production cross section times acceptance times efficiency a x A x e, 
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Table 4 Data and SM background expectation in the signal region for the SR1-SR5 selections. For the SM expectations both 
the statistical and systematic uncertainties are included. In each signal region, the individual uncertainties for the different 
background processes can be correlated, and do not necessarily add in quadrature to the total background uncertainty. 


Signal Region 

SRI 

SR2 

SR3 

SR4 

SR5 

Observed events 

364378 

123228 

44715 

18020 

7988 

SM expectation 

372100 ± 9900 

126000 ± 2900 

45300 ± 1100 

18000 ± 500 

8300 ± 300 

Z(—> vv) 

217800 ± 3900 

80000 ± 1700 

30000 ± 800 

12800 ± 410 

6000 ± 240 

W(—t tv) 

79300 ± 3300 

24000 ± 1200 

7700 ± 500 

2800 ± 200 

1200 ± 110 

W{-+ ev) 

23500 ± 1700 

7100 ± 560 

2400 ± 200 

880 ± 80 

370 ± 40 

W(->- fiv) 

28300 ± 1600 

8200 ± 500 

2500 ± 200 

850 ± 80 

330 ± 40 


530 ± 220 

97 ± 42 

19 ± 8 

7 ± 3 

4 ± 2 

Zj ^*(—> ) 

780 ± 320 

190 ± 80 

45 ± 19 

14 ± 6 

5 ± 2 

tt, single top 

6900 ± 1400 

2300 ± 500 

700 ± 160 

200 ± 70 

80 ± 40 

Dibosons 

8000 ± 1700 

3500 ± 800 

1500 ± 400 

690 ± 200 

350 ± 120 

Multijets 

6500 ± 6500 

800 ± 800 

200 ± 200 

44 ± 44 

15 ± 15 


Table 5 Data and SM background expectation in the signal region for the SR6-SR9 selections. For the SM expectations both 
the statistical and systematic uncertainties are included. In each signal region, the individual uncertainties for the different 
background processes can be correlated, and do not necessarily add in quadrature to the total background uncertainty. 


Signal Region 

SR6 

SR7 

SR8 

SR9 

Observed events 

3813 

1028 

318 

126 

SM expectation 

4000 ± 160 

1030 ± 60 

310 ± 30 

97 ± 14 

Z(—> vv) 

3000 ± 150 

740 ± 60 

240 ± 30 

71 ± 13 

W(— y tv) 

540 ± 60 

130 ± 20 

34 ± 8 

11 ± 3 

W(-»- ev) 

170 ± 20 

43 ± 7 

9 ± 3 

3 ± 1 

W{-> fiv) 

140 ± 20 

35 ± 6 

10 ± 2 

2 ± 1 

W(-+ a*+/0 

3 ± 1 

2 ± 1 

1 ± 1 

1 ± 1 


2 ± 1 

0 ± 0 

0 ± 0 

0 ± 0 

tt, single top 

30 ± 20 

7 ± 7 

1 ± 1 

0 ± 0 

Dibosons 

183 ± 70 

65 ± 35 

23 ± 16 

8 ± 7 

Multijets 

6 ± 6 

1 ± 1 

0 ± 0 

0 ± 0 


Table 6 Observed and expected 90% CL and 95% CL upper limits on the product of cross section, acceptance and 
efficiency, a x A x e, for the SR1-SR9 selections. 


Upper limits on a x A x e [fb] 

Signal Region 

90% CL Observed (Expected) 

95% CL Observed (Expected) 

SRI 

599 (788) 

726 (935) 

SR2 

158 (229) 

194 (271) 

SR3 

74 (89) 

90 (106) 

SR4 

38 (43) 

45 (51) 

SR5 

17 (24) 

21 (29) 

SR6 

10 (14) 

12 (17) 

SR7 

6.0 (6.0) 

7.2 (7.2) 

SR8 

3.2 (3.0) 

3.8 (3.6) 

SR9 

2.9 (1.5) 

3.4 (1.8) 


using the CL S modified frequentist approach [124] and considering the systematic uncertainties on the SM 
backgrounds and the uncertainty on the quoted integrated luminosity. The results are presented in Table 6. 
Values of a x A x e above 599 fb-2.9 fb (726 fb-3.4 fb) are excluded at 90% CL (95% CL) for SR1-SR9 
selections, respectively. Typical event selection efficiencies varying from 88% for SRI and 83% for SR3 to 
82% for SR7 and 81% for SR9 are found in simulated Z(—> vv)+]ets background processes. 
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Fig. 6 Measured distributions of (a) the jet multiplicity, (b) i?™ lss , (c) leading jet fy \-, and (d) the leading jet fy\- to i?™ lss 
ratio for the SRI selection compared to the SM expectations. The Z (—t i/p)+jets contribution is shown as constrained by 
the W(—t //v)+jets control sample. Where appropriate, the last bin of the distribution includes overflows. For illustration 
purposes, the distribution of different ADD, WIMP and GMSB scenarios are included. The error bands in the ratios shown 
in lower panels include both the statistical and systematic uncertainties on the background expectations. 


8.1 Large extra spatial dimensions 

The results are translated into limits on the parameters of the ADD model. The typical Axe of the selection 
criteria vary, as the number of extra dimensions n increases from n = 2 to n = 6, between 23% and 33% for 
SRI and between 0.3% and 1.4% for SR9, and are approximately independent of Mjj. 

The experimental uncertainties related to the jet and E™ lss scales and resolutions introduce, when 
combined, uncertainties in the signal yields which vary between 2% and 0.7% for SRI and between 8% 
and 5% for SR9, with increasing n. The uncertainties on the proton beam energy result in uncertainties on 
the signal cross sections which vary between 2% and 5% with increasing n, and uncertainties on the signal 
acceptance of about 1% for SRI and 3%-4% for SR9. The uncertainties related to the modelling of the 
initial- and final-state gluon radiation translate into uncertainties on the ADD signal acceptance which vary 
with increasing n between 2% and 3% in SRI and between 11% and 21% in SR9. The uncertainties due to 
PDF, affecting both the predicted signal cross section and the signal acceptance, result in uncertainties on 
the signal yields which vary with increasing n between 18% and 30% for SRI and between 35% and 41% for 
SR9. For the SRI selection, the uncertainty on the signal acceptance itself is about 8%-9%, and increases 
to about 30% for the SR9 selection. Similarly, the variations of the renormalization and factorization scales 
introduce a 9% to 30% change in the signal acceptance and a 22% to 40% uncertainty on the signal yields 
with increasing n and E™ 1SS requirements. 

The signal region SR7 provides the most stringent expected limits and is used to obtain the final results. 
Figure 8 shows, for the SR7 selection, the ADD a x A x e as a function of Mp for n = 2, n = 4, and 
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Fig. 7 Measured distributions of the jet multiplicity, leading jet pxi an d the leading jet px to _B™ 1SS ratio for (a) SR7 
and (b) SR9 selections compared to the SM expectations. The Z(—¥ fpj+jets contribution is shown as constrained by 
the W(—t //v)+jets control sample. Where appropriate, the last bin of the distribution includes overflows. For illustration 
purposes, the distribution of different ADD, WIMP and GMSB scenarios are included. The error bands in the ratios shown 
in lower panels include both the statistical and systematic uncertainties on the background expectations. 


n = 6, calculated at LO. For comparison, the model-independent 95% CL limit is shown. Expected and 
observed 95% CL lower limits are set on the value of Mp as a function of the number of extra dimensions 
considered in the ADD model. The CL S approach is used, including statistical and systematic uncertainties. 
For the latter, the uncertainties on the signal acceptance times efficiency, the background expectations, and 
the luminosity are considered, and correlations between systematic uncertainties on signal and background 
expectations are taken into account. In addition, observed limits are computed taking into account the ±lcr 
LO theoretical uncertainty. Values of Mp> below 5.25 TeV (n = 2), 4.11 TeV (n = 3), 3.57 TeV (n = 4), 
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3.27 TeV ( n = 5), and 3.06 TeV (n = 6) are excluded at 95% CL, which extend significantly the exclusion 
from previous results using 7 TeV data [12]. The observed limits decrease by about 6%—8% after considering 
the —la uncertainty from PDF and scale variations in the ADD theoretical predictions (see Table 7 and 
Fig. 9). 

As discussed in Ref. [12], the analysis partially probes the phase-space region with s > Mp, where VI 
is the centre-of-mass energy of the hard interaction. This challenges the validity of model implementation 
and the lower bounds on Mjj, as they depend on the unknown ultraviolet behaviour of the effective theory. 
For the SR7 selection, the fraction of signal events with s > Mp is negligible for n = 2, but increases with 
increasing n from 1% for n = 3 and 6% for n = 4, to about 17% for n = 5 and 42% for n = 6. The observed 
95% CL limits are recomputed after suppressing, with a weighting factor M^/s 2 , the signal events with 
s > Mp, here referred to as damping. This results in a decrease of the quoted 95% CL on Mjj which is 
negligible for n = 2 and about 3% for n = 6 (see Fig. 9). 



Fig. 8 The predicted ADD product of cross section, acceptance and efficiency, a X A X e, for the SR7 selection as a 
function of the fundamental Planck scale in 4 + n dimensions, M />, for n = 2, n = 4, and n = 6, where bands represent the 
uncertainty on the theory. For comparison, the model-independent observed (solid line) and expected (dashed line) 95% CL 
limits on er x A X e are shown. The shaded areas around the expected limit indicate the expected ±1<t and ±2rr ranges of 
limits in the absence of a signal. 


Table 7 The 95% CL observed and expected limits on the fundamental Planck scale in 4 + n dimensions, Mp, as a function 
of the number of extra dimensions n for the SR7 selection and considering LO signal cross sections. The impact of the ± 1 rr 
theoretical uncertainty on the observed limits and the expected ±1 a range of limits in the absence of a signal are also given. 
Finally, the 95% CL observed limits after damping of signal cross section for .s > (see body of the text) are quoted 
between parentheses. 


95% CL limits on Mjj [TeV] 


n extra 

dimensions 

+la (theory) 

95% CL observed limit 

Nominal (Nominal after damping) 

— 1<t( theory) 

95% CL expected limit 
+lcr Nominal — lcr 

2 

+0.31 

5.25 (5.25) 

-0.38 

-0.59 

5.25 

+0.58 

3 

+0.25 

4.11 (4.11) 

-0.33 

-0.38 

4.11 

+0.36 

4 

+0.20 

3.57 (3.56) 

-0.29 

-0.26 

3.57 

+0.25 

5 

+0.17 

3.27 (3.24) 

-0.25 

-0.23 

3.27 

+0.21 

6 

+0.13 

3.06 (2.96) 

-0.19 

-0.20 

3.06 

+0.18 
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Fig. 9 Observed and expected 95% CL limit on the fundamental Planck scale in 4 + n dimensions, Mq, as a function of 
the number of extra dimensions. In the figure the two results overlap. The shaded areas around the expected limit indicate 
the expected Tier and ±2cr ranges of limits in the absence of a signal. Finally, the thin dashed line shows the 95% CL 
observed limits after the suppression of the events with s > M (damping) is applied, as described in the body of the text. 
The results from this analysis are compared to previous results from ATLAS at 7 TeV [12] without any damping applied. 


8.2 Weakly interacting massive particles 

In the following, the results are converted into limits on the pair production of WIMPs. As illustrated in 
Fig. 1, this is done both in the EFT framework and in a simplified model where the WIMP pair couples to 
Standard Model quarks via a Z' boson. 

For each EFT operator defined in Table 1, the limits on M* are extracted from those signal regions 
that exhibit the best expected sensitivity: these are SR4 for Cl, SR7 for Dl, D5, D8, and SR9 for C5, D9, 
Dll. These are translated into corresponding 95% CL limits on the suppression scale M* as a function of 
m x . To derive these lower limits on A/*, the same CL S approach as in the case of the ADD LED model 
is used. The uncertainties on the WIMP signal acceptance include: a 3% uncertainty from the uncertainty 
on the beam energy; a 3% uncertainty from the variation of the renormalization and factorization scales 
and a 5% uncertainty from the variation of the parton-shower matching scale; a 1% to 10% uncertainty 
from uncertainties on jet and E™ 1SS energy scale and resolution; and a 5% to 29% uncertainty due to PDF, 
depending on the operator and WIMP mass. 

Similarly, the uncertainties on the signal cross section are: a 2% to 17% (40% to 46%) uncertainty due 
to the variation of the renormalization and factorization scales in Dl, D5 and D9 (C5 and Dll) operators; 
and a 19% to 70% (5% to 36%) uncertainty due to the PDF for C5, Dll and Dl (D5 and D9) operators, 
with increasing WIMP mass. These theoretical cross-section uncertainties are not considered when deriving 
limits and are not displayed in the plots. A 2% to 9% uncertainty on the cross section, due to the beam 
energy uncertainty, is taken into account. 

The M* limits for five of the operators are shown in Fig. 10 down to WIMP masses of 10 GeV, and 
could be extrapolated even to smaller m x values since there is a negligible change in the cross section or 
the kinematic distributions at the LHC for such low-mass WIMPs. The lcr and 2<r error bands around the 
expected limit are due to the acceptance uncertainties (experimental and theoretical). The effect of the 
beam-energy uncertainty on the observed limit is negligible and is not shown. 

Various authors have investigated the kinematic regions in which the effective field theory approach for 
WIMP pair production breaks down [12—45] . The problem is addressed in detail in Appendix A, where the 
region of validity of this approach is probed for various assumptions about the underlying unknown new 
physics. Here, the EFT framework is used as a benchmark to convert the measurement, and in the absence 
of any deviation from the SM backgrounds, to a limit on the pair production of DM (with the caveat of 
not complete validity in the full kinematic phase space). These are the central values of the observed and 
expected limits in Fig. 10. A basic demonstration of the validity issue is also included in the figure. This is 
done by relating the suppression scale JH* to the mass of the new particle mediating the interaction, JH me< j, 
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Fig. 10 Lower limits at 95% CL on the suppression scale M* are shown as a function of the WIMP mass m x for (a) Dl, 
(b) D5, (c) D8, (d) D9, (e) Dll and (f) C5 operators, in each case for the most sensitive SR (SR7 for Dl, D5, D8, SR9 for 
D9, Dll and C5). The expected and observed limits are shown as dashed black and solid blue lines, respectively. The rising 
green lines are the M* values at which WIMPs of the given mass result in the relic density as measured by WMAP [26], 
assuming annihilation in the early universe proceeded exclusively via the given operator. The purple long-dashed line is 
the 95% CL observed limit on M* imposing a validity criterion with a coupling strength of 1, the red dashed thin lines are 
those for the maximum physical coupling strength (see Appendix A for further details). 


and the coupling constants of the interaction, gi by 

-Mmed = /(Si,M*) . 

For such a relation, an assumption has to be made about the interaction structure connecting the initial 
state to the final state via the mediator particle. The simplest interaction structures are assumed in all cases. 
The form of the function / connecting M med and M* depends then on the operator (see Appendix A). For 
a given operator, one possible validity criterion is that the momentum transferred in the hard interaction, 
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Q tr, is below the mediator particle mass: Qtr < M med . According to this criterion, events are omitted where 
the interaction energy scale exceeds the mediator particle mass. This depends on the values adopted for 
the couplings. Two values (one and the maximum possible value for the interaction to remain perturbative) 
are used. After reducing the signal cross section to the fraction of remaining events, the mass suppression 
scale M* can be rederived yielding potentially two additional expected truncated limit lines in Fig. 10. The 
truncated limits fulfil the respective validity criteria wherever the lines are drawn in the figure. For D9 
for example, the maximum couplings criterion is fulfilled for all WIMP masses, the coupling equal to one 
criterion is fulfilled for WIMP masses up to 200 GeV. For C5 on the other hand, the validity criterion for 
a coupling value of one is violated over almost the whole WIMP mass range, and a truncated limit line is 
only drawn up to a WIMP mass of 10 GeV. 

Figure 10 also includes thermal relic lines (taken from Ref. [40]) that correspond to a coupling, set by 
M*, of WIMPs to quarks or gluons such that WIMPs have the correct relic abundance as measured by the 
WMAP satellite, in the absence of any interaction other than the one considered. The thermal relic line for 
D8 has a bump feature at the top-quark mass where the annihilation channel to top quarks opens. Under 
the assumption that DM is entirely composed of thermal relics, the limits on M* which are above the value 
required for the thermal relic density exclude the case where DM annihilates exclusively to SM particles 
via the corresponding operator. Should thermal relic WIMPs exist in these regions (above the thermal relic 
line), there would have to be other annihilation channels or annihilation via other operators in order to be 
consistent with the WMAP measurements. 

Another way to avoid the validity issues discussed above is to use a simplified model to explicitly 
parameterize the interaction of quarks or gluons with WIMP pairs via generic interactions with real mediator 
particles. With this approach, the coupling of pairs of Dirac fermion WIMPs to quarks via a vector mediator 
particle (such as a Z' boson, corresponding to the operator D5) of a given mass and width (M med and F, 
respectively) is probed. Given the cross-section limit and using simulations at fixed values of M med and F, 
the product of the coupling constants of the Z' boson to quarks and WIMPs, ^/g 9 g x , can be constrained. 
This constraint corresponds to one value in the M*-M med plane as shown in Fig. 11(a), since the mass 
suppression scale can be calculated exactly in this model, M* = M med / % /g7gx. The figure demonstrates 
how, for a given mediator particle mass and two values of the width T, the real value of the mass suppression 
scale would compare to the M* value derived assuming a contact interaction (shown as dashed lines in the 
figure). This contact interaction regime is reached for M med values larger than 5 TeV in the figure. In the 
intermediate range (700 GeV < M med < 5 TeV), the mediator would be produced resonantly and the actual 
M* value is higher than in the contact interaction regime. In this case the contact interaction limits would 
be pessimistic: they would underestimate the actual values. Finally, the small mediator mass regime below 
700 GeV has very small M* limits because the WIMP would be heavier than the mediator, and WIMP pair 
production via this mediator would thus be kinematically suppressed. In this region, the contact interaction 
limits would be optimistic and overestimate the actual M* values. 

In Fig. 11(b) the observed 95% CL upper limits on the product of couplings of the simplified model 
vertex are shown in the plane of mediator and WIMP mass (M med versus m x ). Within this model, the 
regions to the left of the relic density line lead to values of the relic density larger than measured and are 
excluded. 

In the effective operator approach, the bounds on M* for a given m x (see Fig. 10) can be converted to 
bounds on WIMP-nucleon scattering cross sections, which are probed by direct DM detection experiments. 
These bounds describe scattering of WIMPs from nucleons at a very low momentum transfer of the order of 
a keV. Depending on the type of interaction, contributions to spin-dependent or spin-independent WIMP- 
nucleon interactions are expected. As in Ref. [12], the limits are converted here to bounds on the WIMP- 
nucleon scattering cross sections and the results are displayed in Fig. 12. Under the assumptions made in 
the EFT approach, the ATLAS DM limits are particularly relevant in the low DM mass region, and remain 
important over the full m x range covered. The spin-dependent limits in Fig. 12 are based on D8 and D9, 
where for D8 the M* limits are calculated using the D5 acceptances (as they are identical) together with 
D8 production cross sections. Both the D8 and D9 cross-section limits are significantly stronger than those 
from direct-detection experiments. 

The DM limits are shown as upper limits on the WIMP annihilation rate, calculated using the same 
approach as in Ref. [12], in the bottom panel of Fig. 12. The operators describing the vector and axial-vector 
annihilations of WIMPs to the four light-quark flavours are shown in this plot. For comparison, limits on the 
annihilation to uu and qq from galactic high-energy gamma-ray observations by the Fermi-LAT [125] and 
H.E.S.S. [126] telescopes are also shown. The gamma-ray limits are for Majorana fermions and are therefore 
scaled up by a factor of two for comparison with the ATLAS limits for Dirac fermions (see Ref. [12] and 
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(a) (b) 



10' 1 


Fig. 11 (a) Observed 95% CL limits on the suppression scale M* as a function of the mediator mass M rnef j, assuming a 
Z'-like boson in a simplified model and a DM mass of 50 GeV and 400 GeV. The width of the mediator is varied between 
Mmed/3 and M me( The corresponding limits from EFT models are shown as dashed lines; contour lines indicating a 
range of values of the product of the coupling constants (^/gq g x ) are also shown, (b) Observed 95% CL upper limits on the 
product of couplings of the simplified model vertex in the plane of mediator and WIMP mass (M me< j versus m x ). Values 
leading to the correct relic abundance [26] are shown by the black solid line. 


references therein for further discussions and explanations). The annihilation rate that corresponds to the 
thermal relic density measured by WMAP [26] and PLANCK [25] satellites is also shown for comparison in 
the figure. 

Finally, Fig. 12 also demonstrates the impact of the EFT validity and the truncation procedure explained 
above on the quoted upper limits for the WIMP-nucleon scattering and WIMP annihilation cross sections. 
The effect depends strongly on the operator and the values for the couplings considered. In general, the 
limits remain valid for WIMP masses up to 0(100) GeV. The variation of the coupling strengths considered 
leads to changes in the quoted cross-section limits of up to one order of magnitude. 


8.3 Associated production of a light gravitino and a squark or gluino 

The results are also expressed in terms of 95% CL limits on the cross section for the associated production 
of a gravitino and a gluino or a squark. As already discussed, a SUSY simplified model is used in which the 
gluino and squark decays lead to a gravitino and a gluon or a quark, respectively, producing a monojet-like 
signature in the final state. Squark and gluino masses up to 2.6 TeV are considered. The acceptance and 
efficiency Axe for the SUSY signal depends on the mass of the squark or gluino in the final state and 
also on the relation between squark and gluino masses. As an example, in the case of squarks and gluinos 
degenerate in mass (m§ = m-q), the signal Axe for the SR7 (SR9) selection criteria is in the range 25%-45% 
(10%-35%) for squark and gluino masses of about 1-2 TeV. 

The systematic uncertainties on the SUSY signal yields are determined as in the case of the ADD and 
WIMP models. The uncertainties related to the jet and E ™ ISS scales and resolutions introduce uncertainties 
in the signal yields which vary between 2% and 16% for different selections and squark and gluino masses. 
The uncertainties on the proton beam energy introduce uncertainties on the signal yields which vary between 
2% and 6% with increasing squark and gluino masses. The uncertainties related to the modelling of initial- 
and final-state gluon radiation translate into a 10% to 15% uncertainty on the signal yields, depending on 
the selection and the squark and gluino masses. The uncertainties due to PDF result in uncertainties on 
the signal yields which vary between 5% and 60% for squark and gluino masses increasing from 50 GeV 
and 2.6 TeV. Finally, the variations of the renormalization and factorization scales introduce a 15% to 35% 
uncertainty on the signal yields with increasing squark and gluino masses. 

Figure 13 presents, for the SR7 and SR9 selections and in the case of degenerate squarks and gluinos, 
uxAxeasa function of the squark/gluino mass for different gravitino masses. For comparison, the model- 
independent 95% CL limits are shown. For each SUSY point considered in the gravitino-squark/gluino 
mass plane, observed and expected 95% CL limits are computed using the same procedure as in the case 
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Fig. 12 Inferred 90% CL limits on (a) the spin-independent and (b) spin-dependent WIMP-nucleon scattering cross section 
as a function of DM mass m x for different operators (see Sect. 1). Results from direct-detection experiments for the spin- 
independent [127—133] and spin-dependent [134—138] cross section, and the CMS (untruncated) results [14] are shown for 
comparison, (c) The inferred 95% CL limits on the DM annihilation rate as a function of DM mass. The annihilation rate is 
defined as the product of cross section a and relative velocity v, averaged over the DM velocity distribution ({u v)). Results 
from gamma-ray telescopes [125,126] are also shown, along with the thermal relic density annihilation rate [25,26]. 


of the ADD and WIMPs models. This is done separately for the different selections, and the one with the 
most stringent expected limit is adopted as the nominal result. In the region with squark/gluino masses 
below 800 GeV, SR7 provides the best sensitivity while SR9 provides the most stringent expected limits for 
heavier squark/gluino masses. Figure 14 presents the final results. Gravitino masses below 3.5 x 10~ 4 eV, 
3 x 10~ 4 eV, and 2 x 10 -4 eV are excluded at 95% CL for squark/gluino masses of 500 GeV, 1 TeV, and 
1.5 TeV, respectively. The observed limits decrease by about 9%—13% after considering the — la uncertainty 
from PDF and scale variations in the theoretical predictions. These results are significantly better than 
previous results at LEP [54] and the Tevatron [15], and constitute the most stringent bounds on the gravitino 
mass to date. For very high squark/gluino masses, the partial width for the gluino or squark to decay into a 
gravitino and a parton becomes more than 25% of its mass and the narrow-width approximation employed 
is not valid any more. In this case, other decay channels for the gluino and squarks should be considered, 
leading to a different final state. The corresponding region of validity of this approximation is indicated in 
the figure. Finally, limits on the gravitino mass are also computed in the case of non-degenerate squarks and 
gluinos (see Fig. 15). Scenarios with rrig = 4 x rriq, rrig = 2 x rriq, rrig = 1/2 x rriq, and rrig = 1/4 x rriq have 
been considered. In this case, 95% CL lower bounds on the gravitino mass in the range between 1 x 1CF 4 eV 
and 5 x 10~ 4 eV are set depending on the squark and gluino masses. 
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Fig. 13 Cross section times acceptance times efficiency a x Ax e for gravitino+squark/gluino production as a function of 
the squark/gluino mass rriq/g in the case of degenerate squarks and gluinos and different gravitino masses for (a) SR7 and 
(b) SR9, compared with the corresponding model-independent limits. 



Fig. 14 Observed (solid line) and expected (dashed line) 95% CL lower limits on the gravitino mass mg as a function 
of the squark mass mriq for degenerate squark/gluino masses. The corresponding dotted line indicates the impact on the 
observed limit of the — la LO theoretical uncertainty. The shaded bands around the expected limit indicate the expected 
Tier and ±2 <t ranges of limits in the absence of a signal. The region above the red dotted line defines the validity of the 
narrow-width approximation (NWA) for which the decay width is smaller than 25% of the squark/gluino mass. 


8.4 Invisibly decaying Higgs-like boson 

The results are translated into 95% CL limits on the production cross section times the branching ratio for 
a Higgs boson decaying into invisible particles as a function of the boson mass. The SR3 selection provides 
the best sensitivity to the signal and it is used for the final results. The dxeof the selection criteria depends 
on the production mechanism and the boson mass considered. In the case of the gg —>• H process, the Axe 
varies between 0.1% and 0.7% with increasing boson mass from 115 GeV to 300 GeV. It varies between 1% 
and 2% for the VV H production process, and varies between 1% and 12% in the VH case. The gg —> H 
process dominates the signal yield and constitutes more than 52% and 67% of the boson signal for a boson 
mass of 125 GeV and 300 GeV, respectively. 

The uncertainties related to the jet and E™ 1SS scales and resolutions introduce uncertainties in the signal 
yields for the SR3 signal region which vary between 10% and 6% for the gg —» H and VV —> H processes 
as the boson mass increases. Similarly, in the case of VH production processes, these uncertainties vary 
between 8% and 4% with increasing mass. The variations of the renormalization and factorization scales 
introduce a 8% to 6%, 0.2% to 0.8%, and 1% to 3% uncertainty on the boson signal yields for gg H, 
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Fig. 15 Observed (solid line) and expected (dashed line) 95% CL lower limits on the gravitino mass mg as a function of the 
squark mass mq for non-degenerate squark/gluino masses and different squark/gluino mass configurations. The dotted line 
indicates the impact on the observed limit of the —lcr LO theoretical uncertainty. The shaded bands around the expected 
limit indicate the expected dzlcr and ±2tr ranges of limits in the absence of a signal. The region above the red dotted 
line defines the validity of the narrow-width approximation (NWA) for which the decay width is smaller than 25% of the 
squark/gluino mass. 


VV —> H , and VH processes, respectively, as the mass increases. The uncertainties due to PDF result in 
uncertainties on the signal yields which vary between 7% and 8%, 2% and 4%, and 2% and 4% for gg —> H, 
VV — > H , and VH processes, respectively. The uncertainty on the parton shower modelling results in a 7% 
uncertainty in the signal yields for the different channels. 

Figure 16 shows the observed and expected 95% CL limits on the cross section times branching ratio 
a x BR (H —> invisible) as a function of the boson mass, for masses in the range between 115 GeV and 
300 GeV. Values for a x BR(H —> invisible) above 44 pb for ra# = 115 GeV and 10 pb for m # = 300 GeV 
are excluded. This is compared with the expectation for a Higgs boson with BR (H —> invisible) = 1. For a 
mass of 125 GeV, values for a x BR (H —> invisible) 1.59 times larger than the SM predictions are excluded 
at 95% CL, with an expected sensitivity of 1.91 times the SM predictions. This indicates that, for a mass 
of 125 GeV, this result is less sensitive than that in Ref. [58] using ZH(Z —> l + l~) final states, and it does 
not yet have the sensitivity to probe the SM Higgs boson couplings to invisible particles. Nevertheless, for 
a Higgs boson mass above 200 GeV this analysis gives comparable results. 


9 Conclusions 

In summary, results are reported from a search for new phenomena in events with an energetic jet and large 
missing transverse momentum in proton-proton collisions at yfs = 8 TeV at the LHC, based on ATLAS 
data corresponding to an integrated luminosity of 20.3 fb” 1 . The measurements are in agreement with the 
SM expectations. The results are translated into model-independent 90% and 95% confidence-level upper 
limits on a x A x e in the range 599-2.9 fb and 726-3.4 fb, respectively, depending on the selection criteria 
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Fig. 16 The observed (solid line) and expected (dashed line) 95% CL upper limit on a X BR( H —» invisible) as a function 
of the boson mass mj.j. The shaded areas around the expected limit indicate the expected ± I a and ±2 a ranges of limits in 
the absence of a signal. The expectation for a Higgs boson with BR (H invisible) = 1, rj H , is also shown. 


considered. The results are presented in terms of limits on the fundamental Planck scale, Mjj, versus the 
number of extra spatial dimensions in the ADD LED model, upper limits on the spin-independent and 
spin-dependent contributions to the WIMP-nucleon elastic cross section as a function of the WIMP mass, 
and upper limits on the production of very light gravitinos in gauge-mediated supersymmetry. In addition, 
the results are interpreted in terms of the production of an invisibly decaying Higgs boson for which the 
analysis shows a limited sensitivity. 
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Appendix A: On the validity of the effective field theory used to describe dark-matter pair 
production 

Appendix A.l: Introduction 

The effective field theories (EFTs) used here are based on the assumption that a new mediator particle 
couples Standard Model particles to pairs of DM particles and that the mediator particle mass is considerably 
larger than the energy scale of the interaction. In such a case the mediator cannot be produced directly 
in LHC collisions and can be integrated out with an EFT formalism. This heavy-mediator assumption is 
indeed justifiable in direct detection WIMP scattering experiments due to the very low momentum exchange 
typically of order keV in the scattering interactions. This assumption is not always correct at the LHC, 
where the momentum transfer reaches the TeV scale [42-45] . 

A minimal condition for the EFT to be valid is that the momentum transferred in the hard interaction at 
the LHC does not exceed the mediator particle mass, thus ensuring that the mediator cannot be produced 
directly: Qtr < Af me d- To probe this validity, further assumptions have to be made about the actual form 
of the interaction vertex, and thereby about the (unknown) interaction structure itself, connecting quarks 
or gluons to WIMPs. 

The simplest of such assumptions are made below for all the operators used here to derive expressions 
for M med , M*, and the interaction coupling constants, to probe the minimal validity criterion. 


Appendix A.2: Connecting M* to M med 

The simplest interaction structure for the operators D5, D8, and D9 is an s-channel diagram, where the 
mediator particle couples to the initial-state quarks and the final-state WIMPs. This interaction is described 
by three parameters, the mediator mass A/ med , the quark-mediator coupling constant gr q , and the mediator- 
WIMP coupling constant g x . The relation of these parameters to M* is 

Turned \/9<l9xM+ ■ 

The theory is no longer in the perturbative regime if the couplings are outside of the range 0 < y/gqg x < 

47T. 

The simplest s-channel diagram for the operators D1 and Cl involves the exchange of a scalar mediator 
particle where the quark—mediator coupling constant is a Yukawa coupling j/ q . In this case, the mediator 
particle masses can be expressed as: 


D1 

Cl 

m q V<i9x 

Ml ed 

m q Vq\ x v\ 

Mi Ml e d 

Let V\ — CA-^med 

AZ^ed = y/Vq.9x ■ \J/m q 

AI m ed = yqX x £\ ■ M*/m q 

In the above, X x is used for scalar coupling strengths. The vacuum expectation value (VEV) of the 
trilinear scalar vertex is represented by v\. The VEV is then related to the mediator mass scale by v\ = 
Ca Tf mc d> where the common assumption of ss 1 is used. The perturbative range is then 0 < yjy<\g x < 4-7T 
for D1 and 0 < HqX x (\ < (47r) 2 £ A for Cl. 

The operators Dll and C5 describe gluons coupling to WIMPs through a loop diagram, requiring 
different expressions relating M* to M med : 

Dll 

C5 

Qs a s9x 

Qf s QtsXx^X 

4M* 3 Ml ed A s 

4 Ml Ml ed A s 

M med = 

Aimed = 

Let a = 4b _1 

Aimed = ^a.9x M* 

Let a = 4b -1 

Aimed = \J O^xCa AI* 
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Table 8 Relations between the mediator mass M me< j and the suppression scale M* for the simplest interaction vertices 
matching the EFT operators considered here. 


Operator (s) 

Relation between M mec i and M* 

Coupling term range 

D1 

M me d = yjy^9x V M */m q 

o < < 47r 

Cl 

Mmed = VqX x C\ M%/m q 

o < 2 AiWCa < (4tt) 2 Ca 

D5, D8, D9 

^med — y/9q_9x 

0 < V9^9x < 4?r 

Dll 

Mmed = tyaThi M+ 

0 < fy a 9x < v46 it 

C5 

^med — \/ «A x Ca M* 

0 < V aX x(x < 4 \ArCx 


A s = i>M med (b > 1) is another mass suppression scale of the loop connected to the initial-state gluons. 
The coupling terms differ from the other operators, as b > 1 ==> 0 < a < 4. As before, v\ = C,\M med , 

and the assumption of C,\ « 1 is used for C5. The perturbative range for the gluon operators is thus 
0 < < \/16n or 0 < ^/a\ x (x < 4:x/n(\ for Dll and C5 respectively. 

A summary of the different relations between M med and M* for each operator of interest and the 
associated coupling ranges is provided in Table 8. All of the operators have a dependence on a coupling 
term, where the value of these couplings is impossible to know without knowledge of the complete theory. 
Scans over the coupling-parameter space are therefore performed below to quantify valid phase-space regions. 

While these relations were derived for s-channel completions, similar validity arguments can be applied 
to the f-channel as a sum of s-channel operators (see Ref. [44] for further details and caveats). 


Appendix A.3: Regions of Validity 

Given a relation between the mediator mass and the suppression scale, Qti < M med can be evaluated and 
the fraction of events fulfilling this validity criterion can be determined. Two different procedures are then 
followed (which were shown to yield the same results). The nominal procedure is a simple truncation, in 
which the signal cross section is rescaled by the fraction of valid events. With this truncated signal cross 
section, new valid limits on the mass suppression scale are derived, MJ alld . 

The second alternative procedure used to cross-check the simple truncation is an iterative procedure 
that scans through M* until a convergence point is reached. 

1. The starting point is the nominal expected limit on M* assuming 100% validity, named M* xp . M* xp is 
set to M| n before executing step 2 for the first time. 

2. For each step i, obtain the relative fraction of valid events Rm d satisfying Qtr < M“ ed , where M™ ed 
is the mediator mass limit obtained in the previous step (depending on M| n ). 

. r -|l/2(d—4) 

3. Truncate M* following Ref. [43] : M° u = RM med M( n , noting that D1 and Dll are dimension 

d = 7 operators, while D5, D8, D9, Cl, and C5 are dimension d = 6. 

4. Go to step 2, using the current M° ut as the new Ml", repeating until the fraction of valid events at a 
given step R X M reaches 0 or 1. 

5. Calculate the total validity fraction Rm^ = R ^M med an d the truncated limit on the suppression scale 

i 

JV/r lid = [RMl d ] 1 / 2 (d_ 4 ) M: xp . 

The fraction of valid events and the truncated limits on M* can be used to assess the validity of the 
EFT approach. In Figs. 17 and 18, this is shown for D1 and Cl. The majority of the parameter space is 
invalid. The operators D1 and Cl are still valid for regions of parameter space with large coupling values 
and low WIMP masses. 

The validity of the vector, axial-vector, and tensor couplings to quarks via the D5, D8, and D9 operators, 
respectively, are much more justifiable, as shown for D9 in Fig. 19. The operator D9 is valid for the majority 
of parameter space, across couplings and WIMP masses, except for the highest values of m x considered. 
While only D9 is valid for the common canonical choice of p q = g x = 1, the other two operators are valid 
for only slightly larger couplings. 

An assessment of the validity of the gluon EFT operators requires the most assumptions, and has a very 
different coupling range under the assumptions discussed in Appendix A.2. Under these assumptions, Dll 
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Fig. 17 (a) The fraction of valid events and (b) truncated limits for D1 at 95% CL as a function of the WIMP mass m x 
and couplings. The white numbers correspond to the minimum coupling value for which M* alld /Mj xp > 99%. The upper 
perturbative coupling limit for D1 is 47 t. 
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Fig. 18 (a) The fraction of valid events and (b) truncated limits for Cl at 95% CL as a function of the WIMP mass m x 
and couplings. The white numbers correspond to the minimum coupling value for which M* alid /Mj xp > 99%. The upper 
perturbative coupling limit for Cl is (47r) 2 C\, where is taken to be 1. 
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Fig. 19 (a) The fraction of valid events and (b) truncated limits for D9 at 95% CL as a function of the WIMP mass m x 
and couplings. The white numbers correspond to the minimum coupling value for which M* alld /Mj xp > 99%. The upper 
perturbative coupling limit for D9 is 47r. 
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Fig. 20 (a) The fraction of valid events and (b) truncated limits for Dll at 95% CL as a function of the WIMP mass m x 
and couplings. The white numbers correspond to the minimum coupling value for which M™ lld /Mj xp > 99%. The upper 
perturbative coupling limit for Dll is ^16-it. 
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Fig. 21 (a) The fraction of valid events and (b) truncated limits for C5 at 95% CL as a function of the WIMP mass m x 
and couplings. The white numbers correspond to the minimum coupling value for which M™ lld /M* xp > 99%. The upper 
perturbative coupling limit for C5 is , where C,\ is taken to be 1. 


and C5 operators are valid for regions of parameter space with large coupling values and low WIMP masses, 
as shown in Figs. 20 and 21, respectively. 

In general, the validity of the EFT operators is better for low WIMP masses. This is important, as 
collider searches are most competitive with other types of experiments at low m x . Additionally, Figs. 17-21 
show how the truncated limit M* alld quickly approaches the nominal limit M* xp . Some operators have larger 
validity regions than others because the M* limits are larger, and it is thus more likely that Qtr < M me< j. 
Stronger limits therefore remain strong, while weak limits are in fact even further diminished by validity 
considerations. 

Truncated limits are more conservative than the corresponding simplified model used for the completion, 
so long as M me & in the model is greater than or equal to the value used for the truncation. This can be seen 
comparing the D5 operator in Fig. 10 with the corresponding simplified model in Fig. 11. 
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